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1. INTRODUCTION 


In order to extend the analysis capability of the MIRADCOM 
Propulsion Mechanics Function, the TEXGAP-3D static linear elastic 
finite element computer code was obtained in July 1977 to be placed 
on the CDC 6600 computer at Redstone Arsenal. The development of the 
code evolved from a pilot effort in 1973, funded by Thiokol/Huntsville 
Division and more recently (the past two and a half years) has been 
sponsored by the US Air Force Rocket Propulsion Laboratory at Edwards, 
California. A short course on the TEXGAP-3D code, sponsored by Pacifica 
Technology, was held in Los Angeles the latter part of June 1977 at 
which time the code was released. Previous successful employment of 
the TEXGAP-2D computer code (developed by the same research group) by 
personnel of the Propulsion Mechanics Function primarily led to the 
decision to expend the effort necessary to get the TEXGAP-3D code 
operational on the Redstone computer system. 


Once the code was adapted to the Redstone CDC 6600 computer and 
made operational, a task requiring a number of minor modifications to 
the obtained version of the code, a check-out procedure was begun. The 
check-out of the code involved running six classical check problems 
and two additional example problems and comparing the results obtained 
(computer output) with the results of computer runs (of the same prob- 
lems) provided by the developers of TEXGAP-3D. Once these comparisons 
were made it was concluded that the code was operating correctly on the 
Redstone CDC 6600 computer. 


At this point, the utilization of TEXGAP-3D for the purpose for 
which it was acquired was undertaken. The TEXGAP-3D analysis to which 
this report pertains is that of a transition region in a typical solid 
propellant rocket motor. The specific details of the analysis and a 
discussion of the results are given in the subsequent sections of this 
report. 


ll. ROCKET GRAIN CONFIGURATION 


The rocket grain configuration is taken from Reference 1 and 
reproduced, in part, in Figure 1. Since only the transition region is 
to be considered, only that portion of the total geometry is shown in 
Figure 1. The transition is from a four-point star cross-section to a 
circular cylindrical bore cross=section. The region considered actually 
extends somewhat beyond (a distance of 0.463 inches) the transition 
region, as shown in Figure 1. Two slightly different star geometries 
are coneidered, as shown on Section A2 in Figure 1 and labeled I and ITI. 
Due to geometrical symmetry, as noted in Figure 1, and to the symmetry 
of the loading conditions to be considered, it was only necessary to 
model a 45° segment of the grain cross-section in the analysis. 


Wi. «=o FINITE ELEMENT MODELING 


The finite element modeling of the transition region of 
Figure 1 will be discussed in detail and specific reference will be 
made to the TEXGAP-3D user's manual (Reference 2). 


A. Grid Generation 


The three-dimensional finite element grid system for 
the transition region was generated using a BLOCK command (Reference 2) 
to create a three-dimensional block region, such as that shown in 
Figure 2, the edges of which may be defined by a cubic interpolation. 
This permits a high degree of definition of curved boundaries. The 
eight corner points must be input, as per Reference 2, while points 
9 through 32 (see Figure 2) not specified are located automatically by 
linear interpolation between corners. As many blocks as are necessary 
to accurately define a region may be generate. 


B. Element Generation 


The three~dimensional isoparametric 20 node brick element 
shown in Figure 3 was used in this analysis. As many bricks as are 
desired may be generated to fill the previously developed block(s). 
These elements are generated by a BRICK or BRICKH command (Reference 2). 
The BRICKH command employs a reformulated element (Reference 3) which 
takes into account near incompressibility (Poisson's ratio approaching 
0.5). In this analysis the propellent grain elements were generated 
with a BRICKH command and the case elements were generated with a BRICK 
command. For both analysis models I and II, initially 16 propellant 
elements and 8 case elements were generated, as can be noted in Figures 4 
and 5. The result of subsequent "rezoning" can also be noted in 
Figures 4 and 5, and this subject will be discussed a little further 
along in the report. 


Cc. Boundary Conditions 


Boundary conditions may be applied to element faces and/or 
element nodes. This is done through various commands, such as PRESSURE, 
SLOPE, etc. (Reference 2), along with a face number or node number (see 
Figure 3) designation. Various boundary condition commands were util- 
ized in this analysis, and of particular importance were those displace- 
ment boundary conditions, taken from previous two-dimensional analyses 
(see Figure 6), which were applied at each end of the transition region. 
These will be mentioned relative to each of the separate analyses to be 
presented. 


D. Post Processing 


Various post processing commands may be utilized to 
compute stresses, strains and displacements at different points in 
an element. Some of these commands will be demonstrated when each of the 
separate analyses is presented. 


E. Rezoning 


The REZONE command (Reference 2) allows one to subdivide 
a coarse grid region of elements, for which a solution has been obtained, 
into a region of smaller elements, and thereby obtain a more accurate 
solution. A REPOS command (Reference 2) allows one to reposition nodal 
points when a better geometry definition is desired in a rezone analysis. 
The REZONE command was employed in this alalysis and the rezoned regions 
can be observed in Figures 4 and 5. 


The 
for both 


MODEL I 


finite element model I is shown in Figure 4. It was analyzed 
pressure and thermal loading. The TEXGAP-3D computer input 


and output is shown on the following pages. 


The 
computer 


1) 


2) 


3) 


Pressure Loading -— Firing at -53°C (-65°F) 
Propellant 


= 1.31 x 10’Pa 


P (19,000 psi) 
= 0.499 
Pp 


pressure = 1425 psi (9.83 x 10°Pa) 
Corresponding end displacements are shown on computer input 
and referenced to Figure 6. 


Thermal Loading (2 Weeks Storage at -53°C (-65°F)) 


= 4,82 x 10°Pa = 2.07 x 10)°pa 


(700 psi) © (3 x 10 psi) 


= 0.499 ee: 0.3 
0.00089 in/in/°F 


Ty = 60°C (140°F) 
Corresponding end displacements are shown on computer input 
and referenced to Figure 6. 


results of three analyses on Model I are shown on the following 


output sheets: 


pressure loading with two sequential rezoning 
computations, without any repositioning of element 
nodal points; 


pressure loading with two sequential rezoning 
computations, with repositioning of some boundary 
element nodal points; 


thermal loading with two sequential rezoning 
computations, without any repositioning of 
element nodal points. 


) 
| 
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The first analysis, with rezoning near the cylindrical bore end 
of the transition (see Figure 4), yielded a maximum strain of 0.0231 
in element 113 (IJK). The second analysis, with rezoning and reposi- 
tioning near the star end of the transition, yielded a maximum strain 
of 0.0594 in element 339. The third analysis, with rezoning near the 
cylindrical bore end, yielded a maximum strain of 0.200 in element 333. 


The following results consist of the full computer output for the 
first analysis (for element orientation purposes, in the interest of 
the reader), and abbreviated output for the second and third analyses. 
The first analysis is contained on pages 8 thru 94, the second on pages 
95 thru 102, and the third on pages 103 thru 108. 
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| -PIWE IW FORWARD ELIMIWATION = 33.627 


WUNBER OF SECTORS (PRUS) = 1452 
| TIME IM BACKSUOSTITUTION = 2.785 
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| 


30 


vc 


BLOCK OPTION 


2 : = Picante 
ale max MIN MAX 
Xo LOC0EC2t o1000€+21 1 06 8 
Y wel GG0Ee21 o10008¢22 se r) 
Z «- -o 000E*21 o1000E+21 Kk 6 0 
PRINT LEVEL = & 
STRESS POINTS FOR BRICK(DEGEN) ELEMENTS 
POINT Sa s2 s3 
ry -1.000 
e -t. 000 
3 -1, 90 
ry -1.900 
s 1.000 | 
6 1.600 | 
7 1.000 | 
8 1, ¢€€ | 
9 4.060 
STRESS PCINTS FCR PRISH ELEMENTS 
POINT Si $2 $3 S6 
1 1.000 0.000 -1.000 
2 04000 0.000 -1.000 
3 0.000 1,000 -1.000 
6 16660 C.0C® 1.000 9 
5 0.000 0.000 1.000 ; 
6 .00¢ 19000 = §©6 1, 608 
7 2333 -333 «0,000 


STRESS POINTS FOR WEOGE ELEMENTS 
POINT Sa s2 ss 


1 2050 @.000 -1. 900 
2 1006¢ @.660 1.000 
3 1-00C 1.000 -1.0C0 
6 0050 = 8. 000 1.000 
5 1-00€ 6.000 1.060 : 
oJ £eQ0@ 40006 te CER i 
7 $667 o5cO C000 


STRESS POINTS ECR TETRA ELEMENTS 
POINT Sa s2 s3 S& 


0.000 2.000 0.000 
1.000 0.000 0.000 
@.000 1.06 ¢.000 
@.0ce 9.9C8 1.000 
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MODEL II 
The finite element model II is shown in Figure 5. It was analyzed 


for both pressure and thermal loading. The TEXGAP-3D computer input 
and output is shown on the following pages. 


Pressure Loading - Firing at -65°F 


Propellant 


E,* 19,000 psi BE." 3x 10” psi 


Vv = 0, 
P 0.499 


pressure = 1425 psi 
Corresponding end displacements are shown on computer 
input and referenced to Figure 6. 


Thermal Loading (2 Weeks Storage at -65°F) 


E, = 700 psi 
v_ = 0.499 
Pp 


H 


o, = 0.000089 in/in/°F 


T = 140°F 
Corresponding end displacements are shown on computer 
input and referenced to Figure 6. 


The results of two analyses on Model II are shown on the following 
computer output sheets: 


1) pressure loading with two sequential rezoning 
computations, without any repositioning of 
element nodal points; 


2) thermal loading with two sequential rezoning 
computations, without any repositioning of 
element nodal points. 
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The first analysis, with rezoning near the cylindrical bore end 
of the transition (see Figure 5), yielded a maximum strain of 0.0260 | 
in element 113 (IJK). The second analysis, with rezoning near the ; 
cylindrical bore end, yielded a maximum strain of 0.235 in element 113. | 


The following results consist of abbreviated computer output for 


both the first and second analyses. The first analysis is contained on 
pages 111 thru 124; the second analysis is contained on pages 125 thru 


131. 
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Figure 2. Cubic Grid Generation Block 


| The eight corner points must be input (1) thru 8). 
| Points 9 thru 32 not specified are located auto- 
matically by linear interpolation between corners. 


133 


quowatg y2TAgq (TIeAPEND) PPON OZ *E -P4NBTA 


134 


OL 


€ 
438WNN JGON Ea 5 


Y3aSWON 39V4 


Cylindrical 


a 


rad 


f 
TAN: 


i \ 
s \ 
ae 


135 


sition (I) 


4. Finite Element Model of Star to 
Bore Tran 
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Figure 5. Finite Element Model of Star to Cylindrical 
Bore Transition (II) 
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b) THERMAL LOADING CONDITION 


Figure 6. Grain Transition Interface Displacements 
(From Previous Two-Dimensional Analyses) 
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Based upon the analyses reported herein the following conclusions 


CONCLUSIONS 


can be reached: 


1) 


2) 


The TEXGAP-3D code appears to be highly applicable to 
the analysis of stresses and strains in regions of 
complex geometrics. It is definitely a user oriented 
code and would be adaptive as a design tool. 


Along with the attributes described above, the code 
also appears to be quite efficient from an economic 
standpoint, based upon the computer run times of the 
analyses reported. 


es a 


REFERENCES 


1. Engineering Analysis of Selected Advanced Propulsion Concepts - 
III, Technical Report RK-CR-76, US Army Missile Command, Redstone 
Arsenal, Alabama, September 1976. 


2. E. B. Becker and R. S. Dunham, "TEXGAP-3D - A User Oriented Three- 
Dimentional Static Linear Elastic Stress Analysis Program", Vol. I, 
March 1977. 


3. E. B. Becker and R. S. Dunham, "TEXGAP-3D - A User Oriented Three- 
Dimensional Static Linear Elastic Stress Analysis Program", Vol. II, 
March 1977. 


139 


DISTRIBUTION 
No. of Copies 


Chemical Propulsion Mailing List 83 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 12 


Technological Consultants, Inc. 

Attn: Dr. R. A. Heller 3 
P. 0. Box 554 

Blacksburg, Virginia 24060 


Athena Engineering Company 

Attn: Dr. T. L. Cost 1 
1155 Northwood Lake 

Northport, Alabama 35476 


Rohm and Haas Company 

Huntsville Defense Contract Office 

Attn: Dr. H. M. Shuey 1 
Arthur Murray Building 

723-A Arcadia Circle 

Huntsville, Alabama 35801 


Dr. R. M. Hackett 5 
Box 1537 Station B- 

Vanderbult University 

Nashville, Tennessee 37235 


Batelle 

Durham Office 1 
P. 0. Box 8796 

Durham, North Carolina 27707 


DRCPM-MPE, J. W. Allen 1 
-HAE, W. T. Moore 1 
-MDE, G. B. Nicholas 1 
-CFE, C. R. Massey 1 
-PE-E, M. Wright 4 
-VIE, C. J. Arnold 1 
-RSE, A. R. Maykut 1 
-HFE, R. J. Masucci 1 

DRSMI-LP, Mr. Voigt 1 

DRDMI-T, Dr. Kobler 1 
-TK 1 


140 


DRDMI-TKK 
-TKP 
-TBD 
-TI (Record Set) 
(Reference Copy) 


DISTRIBUTION 


No. of Copies 


